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Spectroscopic  temperature  measurements  were  made  at  the 
focal  point  of  imploding  shock  waves  in  the  UTIAS  implosion 
chamber,  which  has  a 20-cm  diameter  hemispherical  cavity. 

The  chamber  was  filled  with  a stoichiometric  H2-O2  gas  mix- 
ture at  different  initial  pressures  (14  ^ 68  atm) . The 
mixture  wus  ignited  at  the  origin  by  an  exploding  wire  gen- 
erating an  outgoing  detonation  wave  which  reflected  at  the 
chamber  wall  as  an  imploding  shock  wave  (gas-runs).  Addi- 
tional experiments  where  an  explosive  shell  of  PETN  was 
plucnd  at  the  hemispherical  wall  were  also  conducted.  The 
shell  was  detonated  by  the  Impact  of  the  reflected  gaseous 
detonation  wave  at  its  surface,  thereby  generating  an  in- 
tense implosion  wave  (explosive-run) . The  temperatures 
were  moasured  at  the  implosion  focus  using  a medium  quartz 
Hilger  spectrograph  with  an  eight  photocell  polychromator 
attachment  over  the  visible  wavelength  range.  The  measured 
radiation  intensity  distributions  were  fitted  to  blackbody 
curves.  The  temperatures  were  10,000  'v  13,000  K for  gas 
runs,  and  15,000  'v  17,000  K for  explosive  runs.  The  con- 
tinuous spectru  from  photographic  film  and  the  measured 
emissivitios,  which  were  very  close  to  unity,  confirmed 
that  the  plasma  was  a blackbody.  Numerical  studios  using 
the  random  choice  method  (RCM)  and  clussical  strong-shock 
theory  wore  used  to  analyse  tho  flows  in  the  entire  range 
of  the  implosion  process.  Ileul-gas  effects  and  radiation 
losses  were  also  considered.  The  results  were  compared 
with  the  experimental  data  and  good  agreement  was  obtained. 

1.  INTRODUCTION 

Tho  UTIAS  explosive-driven  implosion  chamber  has  been  used  in  several  re- 
search areas  since  it  was  conceived  by  Glass  (Ref,  1)  in  the  early  196n,s. 

The  extremely  high  pressures  and  temperatures  generated  at  the  focus  of  implo- 
sions were  utilized  for  driving  projectiles  to  hypervelocities  (Ref.  2),  as  a 
shock-tube  driver  (Ref,  3),  and  to  produce  diamonds  (Ref,  4),  Currently,  ex- 
perimental studies  were  done  on  deuterium-deuterium  fusion  reactions  and 
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measurements  of  the  resulting  neutrons  and  gamma  rays  (Ref,  S) . 

The  c<  editions  produced  by  imploding  spherical  shock  waves  were  first  ob- 
tained by  Gudorly  (Ref.  6).  His  self-similar  solutions,  however,  deal  only 
with  very  strong  shock  waves.  They  do  not  allow  for  any  characteristic 
lengths,  such  as  the  chamber  radius  in  the  present  problem.  Consequently,  in 
the  present  study,  the  conditions  in  the  chamber  wore  calculated  by  combining 
the  numerical  results  of  the  RCM  and  solutions  using  strong-shock  theory. 

The  principle  of  operation  of  the  UTIAS  20-cm  dia  implosion  chamber  can  be 
understood  by  referring  to  Fig,  1,  The  hemispherical  chamber  is  filled  with  a 
2H2+O2  mixture  at  high  pressure  (14  'v  68  atm  for  pure  gas  runs  and  27  atm  for 
explosive  runs).  The  gas  is  ignited  by  exploding  a fine  nickel  wire  (0,13  mm 
dia  x 1 mm  long)  at  the  origin  of  the  hemisphere  creating  an  outgoing  detona- 
tion wave  which  is  reflected  at  the  hemispherical-chamber  wall  and  then  con- 
verges on  the  origin  as  an  imploding  shock  wave.  For  an  explosive  run,  a PETN 
explosive  shell  bonded  against  the  wall  of  a copper-carrier  liner  is  ignit.ed 
upon  reflection  of  the  detonation  wave  and  further  reinforces  the  imploding 
shock  wave,  The  imploding  wave  is  again  reflected  at  the  origin  leaving  behind 
an  extreme  high-pressuro- temperature  region.  The  objective  of  the  present  work 
was  to  measure  the  peak  temperature  of  this  region  and  compare  the  data  with 
analysis. 
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Pig.  1 Schematic  of  implosion  chamber 

wave  dynamics. 

(a)  Outgoing  detonation  wave  in  2H2+O2 

mixture. 

(b)  Detonation  wave  reflects  as  a shock 
wave  (or  initiates  the  hemispherical 
PETN  liner  when  used) , 

(c)  Implosion  phase  onto  the  origin  when 
tho  converging  shock  becomes  very 
strong. 

(d)  a., ..’.Dion  reflects  as  an  exploding 
shock  wavo  leaving  behind  a high- 
pressure-tomperature  region. 


In  previous  experiments,  time-integrated  (Rof.  7)  as  well  as  time-resolvod 
colour-temperatures  of  the  implosion  focus  made  with  large  observing  areas 
(3  x S mm)  wore  found  to  bo  surprisingly  low  (4000  to  5000  K)  (Ref.  8).  In 
tho  present  work,  a much  smaller  observing  area  (0.55  mm  dia)  was  used  and 
timo-resolvod  simultaneous  measurements  were  made  at  eight  wavelengths 
covering  the  visible  region.  Our  temporature  measurements  are  more  precise 
(Ref.  u ) and  are  in  better  agreement  with  expocted  values.  Also,  the 
emissivity  of  the  radiating  plasma  was  evaluated  for  the  first  time. 
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2.  EXPERIMENTAL  ARRANGEMENTS 


The  UTIAS  implosion  chamber  consists  of  two  massive  front  and  rear  plates 
held  together  by  32  large  bolts  (Ref.  9).  A 20-cm  dia  hemispherical  cavity  is 
machined  in  the  rear  plate  supporting  a hemispherical  copper  liner.  An  explo- 
sive PETN  packago  is  manufactured  on  the  inner  surface  of  the  copper  liner  for 
explosive  runs.  The  barrel,  .vhich  holds  the  Plexiglas  observing  window,  the 
gas  inlet  and  the  exploding  wire  are  installed  in  the  front  plate.  For  explo- 
sive runs,  a 20-degree  conical  liner  plate  is  added  to  the  front  plate  to 
guard  the  chamber  from  any  serious  mechanical  damage  due  to  off-focus  implo- 
sions (Ref,  10). 

A 0.13-mm  dia  x 1-mm  long  nickel  wire  is  exploded  at  the  geometrical  centre 
of  the  chamber  by  discharging  a 1-pF  capacitor  charged  to  21  KV.  All  spectro- 
scopic observations  were  made  with  a llilgor  medium-quartz  spectrograph.  For 
photographic  work  Kodak  Plus-X  pan  film  was  used.  For  time-resolved  photo- 
electric recording  a photomultiplier  attachment  enabled  simultaneous  measure- 
ments at  eight  wavelengths  (3900,  41-8,  4391,  4710,  5123,  5629,  6328  and 
7525  X).  The  time-response  of  tlio  system  was  about  0.1  usee.  The  calibration 
of  the  photomultipliers  was  done  in  the  conventional  manner  using  a calibrated 
tungsten  ribbon  strip  lamp. 

For  gas  runs,  a 2H2+02  gas-mixturo  was  filled  at  different  pressures  (13.6, 
27.2,  40.9,  54.5  and  68,1  atm).  Thickness  of  1.6  mm  and  3.2  mm  PETN  explosive 
sholls  were  used  with  a 2H2+O2  gas-mixturo  of  27.2  atm  initial  filling  pressure 
for  explosive  runs.  These  oxplosivo  shells  woighed  approximately  40g  and  80g, 
respectively. 


3.  EXPERIMENTAL  RESULTS 

Pictures  of  the  radiation  spectrum  of  the  first  implosion  pulse  wore  taken 
to  investigate  the  characteristics  of  tho  radiation  with  a rotating-disk 
shutter.  Uotails  of  tho  operation  of  the  shutter  can  be  found  in 'Refs.  7 and 
11.  It  was  found  that  tho  spoctrum  was  completely  continuous  suggesting  that 
the  temperature  of  tho  radiating  gas  at  the  implosion  focus  can  be  determined 
by  finding  a blackbody  curve  of  a certain  temperature  which  fits  best  to  tho 
experimental  data.  Tho  temperature  measurements  wore  carried  out  photoeloc- 
tritally  with  tho  experimental  sot-up  shown  in  Fig,  2.  A brass  mask  with  a 
small  hole  (0.55  mm  din)  at  the  centre  was  placed  on  the  window  surfaco.  The 
observing  area  was  determined  by  tho  hole  size  and  the  width  of  the  entrance 
slit  of  tho  spectrograph. 


Iropiouor  Chombef 


Fig.  2 Schematic  diagram  of  experi 
mental  facility  to  study  spherical 
implosions. 
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A typical  oscillogram  of  a time-resolved  measurement  for  a gas  run  is  shown 
in  Fig.  3.  The  recordings  are  strongly  disturbed  by  the  electromagnetic  pick- 
up from  the  discharge  system  of  the  exploding  wire.  The  interference  disap- 
pears in  about  30  nsec  and  the  radiation  intensity  stays  at  a quite  constant 
value  until  the  first  implosion  wave  comes  back  to  the  origin.  The  implosion 
time  (the  time  between  the  ignition  and  the  arrival  of  the  implosion  shock  wave 
at  the  centre)  was  approximately  75  nsec  for  gas  runs  and  59  nsec  and  56  nsec 
for  explosive  runs  of  1.6  and  3,2  mm  PETN  layers,  respectively.  In  gas  runs 
some  subsequent  implosion  pulses  were  observed  although  the  radiation  intensi- 
ties were  much  lower  than  the  first  implosion  pulse.  The  window  broke  about 
500  'v  700  nsec  after  ignition  for  gas  runs  and  about  18  nsec  after  the  first 
implosion  pulse  for  explosive  runs. 

Fig,  3.  Typical  output 
signal  from  photomulti- 
plier for  a gas  run. 
po  * 27.2  atm,  X - 5123A. 

At  Noise  from  discharging 
system  ♦ radiation 
from  exploding  wire. 

B:  Implosion  pulse. 


The  radiation  intensities  at  eight  different  wavelengths  which  cover  the 
whole  visible  radiation  range  were  measured  for  each  run  to  eliminate  the  run- 
to-run  fluctuation  due  to  the  variation  in  the  degroo  of  focussing. 


The  radiation  intensity  of  a graybody  at.  a given  temporaturo  T and  wave- 
length Xi  is  expressed  as 


I\i  - B-BxiCT) 


CD 


whore  Bx^CT)  is  the  Planck  law. 
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with  Ci  « 1,191  x 10"  erg.cm  «soc  and  C2  ■ 1.4388  cm  K;  Ii  is  tho  omissivity 
of  the  radluting  gas.  Under  the  condition  expCC^/XT)  >>  1,  B^CT)  can  be  re- 
placed by  tlie  I'ien  approximation, 


Cl  / “c2  \ 

»Xi(T)  *'  — - exp  f j (3) 
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Substituting  Eq.  (3)  into  Eq,  (1)  results  in 
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(4) 


Obtaining  tho  intensities  IX^  (i  ■ 1,  2,  8)  experimentally,  the  tempera- 

ture can  be  determined  from  the  slope  of  the  line  which  goes  through  the  data 
points  of  a plot  of  the  values  of  the  left-hand  side  of  Eq.  (4)  against  tho 
wavo  number  (l/\i).  Tho  emissivity  can  also  be  evaluated  from  tho  values  at 
the  intersection  of  the  line  with  the  vortical  axis.  However,  after  evaluating 
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some  experimental  data,  it  wa3  found  that  the  tomperature  range  obtained  at  the 
implosion  focus  was  above  10,000  K and  the  accuracy  of  Wien's  approximation  was 
poor.  Therefore,  in  practice,  tho  temperature  and  emissivity  were  varied  to 
minimize  tho  sum  of  the  squares  of  deviations  in  tho  experimentally  obtained 
intensities  from  tho  intensities  givon  by  Eq.  (23 , Typical  results  of  a gas 
run  and  an  explosive  run  are  shown  in  Pig.  4.  The  results  are  presented  taking 
advantage  of  Eq,  (4),  namely,  the  solid  lines  are  drawn  in  such  a way  that 
their  slopes  and  their  intersections  with  the  ordinate  represent  the  tempera- 
ture and  emissivity  obtained  by  using  Planck's  law  and  the  data  points  are 
correspondingly  corrected  by  a factor 

C2 

Px±  (T3  “ 1 - exp  (5) 

which  is  simply  tho  ratio  of  Planck's  law  to  Wien's  approximation. 


Fig.  4.  Experimental  rosults  for  a gas  run  and  an  explosive  run. 

It  was  found  that  the  emissivity  of  tho  radiating  gas  was  very  close  to 
unity  except  for  the  low  initial  pressure  gas  runs  (p0  ■ 13.6  and  27.2  atm) 
whore  tho  emissivitios  have  a range  of  0,7  ^ 0.9. 

4.  ANALYTICAL  PREDICTIONS 

In  the  present  study,  the  RCM  was  employed  to  calculate  the  imploding 
stugo.  Details  of  tho  RCM  aro  available  in  Refs.  12  to  IS.  In  the  RCM,  the 
wavos  are  propagated  statistically  and  the  discontinuities  such  as  shocks  and 
contuct  surfaces  appear  as  perfectly  sharp  fronts  without  any  smearing,  unlike 
finite-difference  schemes  using  artificial  viscosity.  Self-similar  solutions 
for  constant  speed  spherically  expanding  Chapman-Jouguet  (C-J)  detonation  waves 
in  the  detonation  stage  wore  calculated  to  give  the  initial  conditions  for  the 
imploding  stage  (Ref.  16).  The  details  of  the  numerical  results  from  -the  RCM 
and  strong-shock  theory  for  gas  runs  and  explosive  runs  aro  given  in  full  in 
Ref.  11. 

5.  COMPARISON  OP  EXPERIMENTAL  RESULTS  AND  THEORY 

figure  S compares  tho  experimental  results  and  the  expected  average  tom- 
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peratures  over  a 0,55-mm  dia  circular  area  at  the  origin.  The  calculated  tem- 
perature depends  on  the  position  of  the  imploding  shock  wave.  When  the  wave 
radius  is  larger  than  the  radius  of  the  obsorving  area  only  tho  wave  front  can 
be  seen  through  the  window.  When  tho  imploding  shock  wave  comes  insido  the 
observing  area,  not  only  the  wave  front  but  also  the  region  behind  the  implodr 
ing  shock  wave  is  seen  and  the  averaged  temperature  becomes  lower  than  the 
temperature  at  the  wave  front.  This  arises  from  the  rapidly  falling  tempera- 
ture behind  the  Implosion  wave.  As  a result,  the  averaged  temperature  has  a 
maximum  when  the  radius  of  the  imploding  shock  wave  is  slightly  smaller  than 
the  radius  of  the  observing  area.  The  imploding  shock  radius  at  which  this 
occurs  is  called  R*,  as  shown  on  Pig.  S.  In  other  words,  there  is  a certain 
maximum  observable  temperaturo  for  a certain  size  of  observing  area  no  matter 
how  small  the  imploding  shock  wave  converges.  This  maximum  temperature  and 
the  temperatures  which  are  expected  to  be  observed  when  the  imploding  waves 
are  at  different  radii  from  the  centre  are  also  shown  in  Pig.  5.  In  principle, 
if  it  wore  possiblo  to  place  tho  obsorving  area  exactly  at  the  focus  of  the 
Implosion  then  higher  temperatures  would  bo  moasurad  with  decreasing  mask  hole 
size.  However,  in  practice  the  size  of  0,55  mm  dia  was  found  to  bo  an  optimum 
limit  considering  the  time  resolution  of  tho  present  measuring  system. 


i 
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Pig.  5(a).  Comparison  of  experimental  Fig.  5(b).  Comparison  of  experimental 
rosults  with  analysis  for  gas  results  with  analysis  for  oxplo- 

runs.  sive  runs, 

For  gas  runs  the  dependence  of  measured  temperaturos  on  tho  initial  filling  J 

pressures  agreed  quite  well  with  the  analytical  predictions  although  the  ab-  j 

solute  values  of  the  temporaturos  were  somowhat  (^  20%)  lower  than  tho  expected 
maximum  temperatures.  Tho  lower  values  of  the  measured  tomperaturo  can  bo 
attributed  to  the  fact  that  the  temperature  behind  the  imploding  shock  front 
decreases  rapidly.  When  the  temperature  is  measured  in  a radial  direction  j 

from  a certain  point  ahead  of  the  wave  front,  it  would  bo  lower  than  tho  front  ; 

temperature.  It  would  probably  be  a kind  of  averaged  temperature  of  tho  region  , 

between  the  shock  front  and  a certain  point  behind  the  front  by  about  the  same 
distance  as  tho  radiation  mean  free  path.  It  has  been  calculated  that,  when 
the  imploding  wave  radius  R«  has  converged  to  R*  ('v  0.25  mm,  in  Pig,  5),  the 
temporatures  at  points  5 and  63u  behind  the  shock  front  are,  respectively,  90 
and  60%  of  tho  shock- front  temperaturo. 

For  the  explosive  run3,  tho  measured  temperaturos  were  much  lower  than  tho 
expected  temperatures.  Also  the  increase  in  the  measured  temperatures  for 
incroased  Pt!TN  thickness  was  small  comrured  with  the  analytical  prediction. 
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Considering  the  much  higher  temperatures  in  explosive  runs  at  the  shock  front 
than  in  gas  runs,  the  so-called  "screening  effect"  of  the  preheated  gas  layer 
nil  cad  of  the  shock  front  (Rof.  17)  is  probably  the  reason  for  this  phenomenon. 

Although  the  measured  maximum  temperatures  were  relatively  low  and  corre- 
sponding to  large  analytical  radii  for  the  imploding  waves  in  explosive  runs, 
the  pressures  obtained  In  explosive  runs  are  much  higher  than  in  gas  tuns 
oven  for  the  highest  initial  gas  pressure  (p0  ■ 68.1  atm)  as  the  window  breaks 
much  fustor  in  explosive  than  in  any  gas  runs.  The  final  size  of  the  implo- 
sion focus  which  gives  the  measured  temperature  obtained  analytically  is  about 
1.5  mm  radius.  The  final  pressure  that  would  be  obtained  in  on  exnlosivu  run 
for  this  focal  size  corresponds  to  the  pressure  which  would  be  ol.f  nod  in  a 
gas  run  of  p0  54  atm  when  it  converged  to  about  0,6  mm  radius,  however, 
the  window  breaks  nearly  20-fold  faster  in  the  explosive  runs,  producing  a 
pressure  far  beyond  what  can  be  produced  in  a gas  run  of  54  atm  initial  pres- 
sure. Despite  tho  fact  that  the  obsorved  temperature  remains  at  a low  value, 
owing  to  the  preheating  of  the  gas  ahead  of  the  shock  wave  and  becomes  opaque 
with  a surface  temperature  of  about  17,000  K,  the  actual  implosion  focus 
continues  to  decrease  to  much  smaller  values  (from  the  1,5  mm  radius).  How- 
ever, this  has  no  effect  on  the  measured  temperatures  as  the  preheated  gas  is 
opaque, 

Although  an  appropriate  perfoct-gas  analysis  (y  'v  1.14)  was  used  through- 
out tho  present  work,  the  approach  appears  quite  reasonable.  Flagg  (Ref.  10) 
has  shown  that  the  isentropic  exponent  of  the  gas  scarcely  changes  during  tho 
implosion  and  reflection  stages  (y  *v  1.14)  after  investigating  the  numerical 
results  of  Brode  (Ref.  18)  who  used  an  equation  of  stute  for  a real  gas  for 
u burnt  2ll2+02  mixturo,  Tho  equilibrium  compositions  of  the  2H2+O2  system 
were  calculated  separately  in  tho  present  study  using  the  method  of  minimizing 
Gibbs'  freo  onorgy  (Rof,  19).  It  was  noted  that  the  offocts  of  increasing 
temperature  and  prossure  on  dissociation  ond  ionization  compensate  oach  other 
and  no  significant  change  in  tho  composition  occurred. 

Tho  radiation  loss  has  also  been  estimated  and  it  was  found  that  the  order 
of  tho  loss  compared  to  the  onorgy  flux  into  the  shock-wave  region  due  to  tho 
converging  geometry,  l.o,,  the  inhomogeneous  term  of  the  energy  equation  in 
tho  lossloss  basic  equations  was  nogligibiy  small  (10-Z  'v  10-6),  owing  to  the 
very  high  pressure  and  particle  velocity  {p  'v  1.4  - 5.8  x 104  atm  and  u a,  12-24 
km/soc  ut  R9  « 0.1  mm).  Therefore,  tho  radiation  loss  .'s  not  important  for 
the  present  work. 

In  tho  most  reconc  experiments  (Rof.  5),  whore  deuterium  was  used  instoud 
of  hydrogen,  production  of  neutrons  and  y-rays  from  D-D  reactions  wore  ob- 
sorved, If  this  is  tho  case  thon  tho  final  focus  of  the  implosions  could  be 
very  small,  say  tho  order  of  lOy. 

6.  CONCLUSIONS 


The  continuous  radiation  spoctra  and  the  measured  emissivities  showed  that 
tho  hot  gases  at  an  implosion  focus  from  gas  and  explosive  drivers  radiate  as 
blackbodies.  The  temperatures  at  tho  focal  point  of  an  imploding  shock  wave 
woro  measured  spectroscopically  in  the  visible  region.  For  gas  runs,  the 
effect  of  tho  initial  filling  pressures  on  the  temperature  was  small  and 

10.000  'v  13,000  K woro  moasured  when  the  initial  pressure  was  increased  from 
13.6  to  68.1  atm.  This  dependence  of  the  temperature  agreed  very  well  with 
analytical  results.  For  explosive  runs,  the  observed  temperatures  wore 

15.000  K for  a J,6-mm  and  17,000  K for  a 3.2-mm  shell  of  PKTN  explosive,  Tho 
measured  tomporaturos  in  explosive  runs  were  limited  to  15,000  'v  17,000  K 
owing  to  tho  scroonlng  effect  of  tho  preheated  layer  ahead  of  the  shock  front, 
even  though  tho  imploding  shock  wovo  continued  to  converge.  A high  degree  of 
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convergence  wus  obtained  despite  the  possible  existence  of  bifurcation  of  the 
Imploding  shock  waves  due  to  the  boundary  layer  developed  on  the  flat  surface 
and  somo  of  tho  unavoidable  geometrical  roughness  of  the  ehambor  around  the 
origin  (which  is  not  small,  if  wo  talk  about  a final  focus  in  tho  micron 
range).  It  is  even  more  remarkable  if  neutrons  and  v-rays  can  be  obtained 
(Ref,  S),  Many  additional  investigations  are  still  needed  to  settle  the 
question  of  the  actual  final  focal  size  and  the  physical  conditions  which 
exist  there,  further  details  can  be  found  in  Ref.  11, 
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